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ABSTRACT 

(N 
> 

Recent laboratory experiments on interstellar dust analogues have shown 

o 

' for the observed H2 abundance in cold diffuse and dense clouds. However, 



o 



X 



that H2 formation on dust grain surfaces is efficient in a range of grain temper- 
atures below 20 K. These results indicate that surface processes may account 



high abundances of H2 have also been observed in warmer clouds, including 
photon-dominated regions (PDRs), where grain temperatures may reach 50 
K, making the surface processes extremely inefficient. It was suggested that 



C/3 

' this apparent discrepancy can be resolved by chemisorption sites. However, 



recent experiments indicate that chemisorption processes may not be efficient 
at PDR temperatures. Here we consider the effect of grain porosity on H2 for- 
mation, and analyze it using a rate equation model. It is found that porosity 
extends the efficiency of the recombination process to higher temperatures. 
This is because H atoms that desorb from the internal surfaces of the pores 
may re-adsorb many times and thus stay longer on the surface. However, this 
porosity-driven extension may enable efficient H2 formation in PDRs only 
if porosity also contributes to significant cooling of the grains, compared to 
non-porous grains. 
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1 INTRODUCTION 

H2 is the most abundant molecule in the Universe. It plays a key role in interstellar chem- 
istry and serves as a coolant during gravitational collapse and star formation. Therefore, 
the processes of formation and dissociation of molecular hydrogen strongly affect the evo- 
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lution of interstellar clouds (|Dulev fc William; 



1984; 



Willi ami 



1993). Already in the 1960's 



it was found that H2 cannot form efficiently enough in the gas phase to account for its 
obser ved abundances, and it was suggested that dust grains serve as catalysts for this pro- 
cess (|Gould fc Salpeteiill963h . The formation of molecular hydrogen on dust grain surfaces 
was studied, assuming that H atoms that collide with a grain are physically adsorbed (ph- 
ysisorbed) . These atoms diffuse on the surface and form molecules when they encounter 



each other (|Hollenbach &: Salpete Jl97nl : lHollenbach fc SalDetedll97lUHollenbach et al.lll97lf) . 

The formation rate i?H2 (cm~^s~^), of H2 molecules can be expressed by 

where nn (cm~'^) is the number density of hydrogen atoms in the gas, fn ■ 



I 



(1) 



(cm s~^) is their typical velocity (where rriH = 1.67-10"^'' gram is the mass of an H atom and 
Tg is the gas temperature), n^^i „ (cm~^) is the number density of dust grains and a (cm^) is 
the average grain cross section (jHoUenbach et al.lll97l[ ). Assuming spherical grains, cr = vrr^ 
where r (cm) is the grain radius. The parameter ^ ^ 1 is the sticking coefficient and 
the parameter ^ r/ ^ 1 is the recombination efficiency, namely, the fraction of adsorbed 
atoms which come out in molecular form. 

The total grain mass amounts to about 1% of the hydrogen mass in the cloud. The 



number density of dust grains is n 



grain 



10 where n = nn + 2nH2 ( cm ^) is the total 



densi ty of hydrogen atoms in both atomic and molecular (riHa) forms (jHoUenbach et al 



19711 1. Assuming compact grains with mass density oi Pg = 2 (gram cm ^), these parameters 



lead to a typical grain diameter of about 0.2/im. Assuming gas temperature of 100 K as well 
as ^ ~ 1 and f] — 1, one obtains 

(2) 



R 



Ha 



Rn^in, 



where the rate constant is R 



clouds 



( jiiral 197,4 



in agreement with observations for diffuse 



More recent observations in photon dominated regions, indicate that in 



these regions the rate co efficient should be R ~ 10 (cm^ s i), in order to account for the 



observed Ho abundance (Habart et al 



20031 ). 



Recent laboratory experiments indicate that a model based on physisorbed atoms may 
account for the H2 formation rate in cold interstellar clouds, where the grain temperature 
does not exceed 20 K. However, molecular hydrogen has also been observed in large abun- 
dances in photon-dominated regions (PDRs), which are exposed to high fluxes of untra- violet 



radiation (|Tielens fc Charnlev 



19971 ). In these regions, typical grain temperatures are in the 
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range of 30 - 50 K. In order to explain the formation of H2 in PDRs it was suggested that 
some of the H atoms are chemically adsorbed (chemisorbed), and thus may remain on the 
surface long enough to form molecules even at higher grain temperatures. However, recent 
experimental and theoretical results indicate that such chemisorption processes are likely to 
contribute to H2 formation only at grain temperatures above 450 K. Such high grain tem- 
peratures can be achived by ultra-small gra ins with diameters of the order of 1 nm, whic h 
may exhibit large temperature fluctuations (jPraine &: AndesronI : iGuhathakurta fc Drainel ). 
However, the efficiency of H2 production on such extremely small grains is expected to be 
low. This is due to the very small population of H atoms on each grain, where the recombi- 
nation ra te is dominated by the di screteness and fluctuations in the population of adsorbed 



H atoms (jLipshtat fc Biham 



20051 ). Thus, chemisorption processes may be irrelevant to the 



range of grain temperatures observed in PDRs. 

In this paper we consider the effect of porosity on the rate of H2 formation on interstellar 
dust grains. We present a model in which H atoms can be adsorbed either on the external 
surface of the grains or inside the pores. The model shows that H atoms can be retained in 
pores for longer times, and recombine there to form H2 molecules at higher temperatures 
than on non-porous grains. 

The paper is organized as follows. In Sec. |21we consider the formation of H2 molecules 
on compact grains, which is efficient only within a narrow temperature range below 20 K. 
In Sec. 0] we describe the high H2 abundance in PDRs, which cannot be explained by the 
processes considered in Sec. El In Sec. 0] we review some of the evidence that interstellar 
dust grains are porous. A model for H2 formation on porous grains is presented in Sec. 
It is shown that porosity extends the efficiency range of H2 formation toward higher grain 
temperatures. The results are discussed in Sec. IHl and summarized in Sec. d 



2 H2 FORMATION ON DUST GRAINS WITH COMPACT SHAPE 

During the past decade, the formation of molecular hydrogen on analogues of interstellar 
dust ma terials has been studie d by laboratory expe rin aents. In particular, carbonaceous 



samples 


(Pirronello et al. 


1997a b: Vidali e 


t al. 


1998alb: 


Pirronello et al. 


1999: Zecho et al. 


2002a 


k 


Perrv & Price 


2f 


03 


Guttler et al. 


2004 


: Zecho et al. 2004). silicates such as olivine 


(Pirronello et al. 


1997a 


b; 


Vic 


lah et al.ll998 


J) and a variety of amorphous ice samples ( Manico et al. 


2001; 


Roser et al. 


2001, 


2002 


Hornekaer et al.l2003. 2005 


Perets et alLoO 


5 ) have been stud- 
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ied. The results of some of these experiments were analyzed using rate equation models. The 
relevant parameters, namely the energy barriers for the atomic diffusion and desorption and 
for molecular desorption, have been found. These parameters were used in order to calculate 
the efficiency of molecular recombination under conditions relevant to the interstellar clouds. 
It was found that the recombination of H atoms physisorbed on these surfaces is efficien ; 



within a narrow t emperature range, below 20 K (jKatz et al 



Perets et al 



20051). 



1999; 



Cazaux fc Tielens 



2004; 



Consider an interstellar dust grain exposed to a flux of H atoms. The formation of 
molecular hydrogen on the grain surface is described by 

^ = F -WN - 2AN^, (3) 
at 

where N is the number of H atoms on the grain. The first term in Eq. Q describes the 
effective flux, F (s^^), of H atoms sticking to the surface, given by F = n-^VYia^. The 
Langmuir-Hinshelwood (LH) rejection of atoms deposited on top of already adsorbed hy- 
drogen atoms (or molecules) is neglected here, since the coverage is assumed to be low. The 
second term describes the thermal desorption of H atoms from the surface where 



W = uexp{-Ei/kBT) 



(4) 



is the desorption rate, z/ is the attempt rate [standardly taken as 10^^ (s~^)], Ei is the 
activation energy barrier for desorption and T (K) is the grain temperature. The third term 
in Eq. accounts for the depletion of adsorbed H atoms due to formation of H2 molecules. 
The parameter A = a/ S is the rate at which H atoms scan the entire surface of the grain, 
where 



uexp{-Eo/kBT) 



(5) 



is the hopping rate between adjacent adsorption sites and Eq is the energy barrier for hop- 
ping. For a spherical grain the number of adsorption sites on the surface is 5* = Arcr'^s where 
s (cm~^) is their density. The formation rate, Rgram, of H2 on a single grain is given by 



R 



gram 



AN"^. Here we assume, for simplicity, that all H2 molecules desorb from the surface 



upon formation. Under s teady state conditi ons one can obtain an exact solution for i?. 



in terms of F, W and A fjBiham et al 
2R, 



1998|). The recombination efficiency 



V 



i-gram 



grain 



(6) 



turns out to exhibit strong temperature dependence. Taking into account the Langmuir 
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rejection term, it is found that H9 recombination is efficient at grain temperatures in the 



range Tq ^ T ^ Ti where (jBiham fc Lipshtatll2n02f ) 



To-- 
and 
Ti = 



kBHuS/F) 
2Ei — Eq 



(7) 



(8) 



kBHuS/F)' 

The simulations presented in this paper are based on the parameters obtained experimentally 
for amorphous car bon, namely the activation energies are Eq = 44.0 meV and Ei = 56.7 meV 



Katz et al 




19991). and the density of adsorption sites on the surface is s ^ 5 x 10 (cm ) 



l.l l2001f ). For simplicity, the sticking coefficient is assumed to be ^ = 1. The gas 
phase density is nn = 100 (cm~^) and its temperature is Tg = 500 K, thus vh = 3.24 x 10^ 
(cm s~^). Under these conditions one obtains that for amorphous carbon the range of grain 
temperatures in which H2 formation is efficient is 12 ^ T ^ 16 K. 



3 THE ABUNDANCE OF H2 IN PHOTON DOMINATED REGIONS 



The temperatures of both the gas and dust grains in PDRs are significantly higher than in 
other H I regions. This is due to far ultraviolet radiation which penetrates th e clouds, heats 



the g as through photon adsorption and the dust through photoelectric heating (jHoUenbach fc Tielens 



19971 ). The temperatures of interstellar grains in several PDR environments, suc h as Chamaele on, 



Oph W, S140, IC63, NGC2023 and the Orion Bar, were recently c alculated 



Habari 



2004^ . An analytic expression for the temperatures of large grains (jHoUenbach et al.lll99l[ ) 



was used, and the median temperature of small grains was evaluated. It was found that in 
most of these regions the average temperatures of small grains are in the range of 10 - 30 K 
(though temperature fluctuations may heat them to temperatures higher than 100 K for very 
short durations), while the temperatures of large grains are between 30 - 50 K. Observations 
show that in all these PDRs, the rate constant of H2 formation is 3 • 10^^'' ^ i? ^ 1.5 ■ 10^^^ 
(cm^ s^^)- This is a rather high formation rate, which is comparable with that obtained in 
colder, diffuse clouds. In order to achieve such rate, the recombination efficiency must be 
very high, namely ~ 1. 

The process of II2 formation on compact grains described above is not efficient at these 
high grain temperatures and cannot account for the observed abundance. This is because at 
these temperatures the residence time of H atoms on the grain is too short to enable them to 
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encounter each other and recor nbine. In order to explain th e formation efficiency at higher 
temperatures it was suggested (jHoUenbach fc SalpetenllQTlh that the grain surface contains 
enhanced adsorption sites with "semi-chemical binding" sites where atoms can stick much 
stronger to the surface. As a result they retain on the surface for a longer period, that enables 
more efficient recombination. It was suggested that these enhanced binding sites are chemi- 
cal adsorption (chemisorpt ion) sites. Their binding energies were calcu l ated, mainly for car 



bonaceous dust materials ( Bennet et al 



1971 



Dovesi et al 



1976, 



1981 



Aronowitz &: Chane 



19851 ). However, recent theoretical and experimental results indicate that chemisorption sites 



play an important role only at much higher surface temperatures. Theoretical studies have 
shown that H atoms on graphite need to pene trate an energy barrier of ^0.2 eV in order 



to reach a chem isorption site on basal planes (jJeloaica fc Sidis 



Sha et al 



1999; 



Sha fc Jackson 



2002; 



20021 ). This result was confirmed experimentally when an activation energy barrier 



of 0.18 eV was found experimentally for hydrogen on C(OOOl) (Zecho, private communica- 
tion). The gas temperatures prevalent in PDRs may enable H atoms from the gas phase 
to enter chemisroption sites upon collision with a grain. However, it is unlikely for atoms 
already adsorbed in physisorption sites (and thus equihbrated with the surface) to penetrate 
this barrier and hop into chemisorption sites. 



Chemisorbed H atoms may reside in defect sites, in step edges (for which the energy 
barrier for penetration is negiligible) or in basal plane sites. The recombination of these 
atoms takes place through two possible processes. In the LH process, a desorbing H atom 
picks another H atom from an adjacent site to form a molecule. In the Eley-Rideal (ER) 
mechanism, an atom from the gas phase abstarcs an adsorbed atom to from a molecule. 
However, the ER mechanism is relevant only at high surface coverage, which is not likely 
to exist under interstellar conditions. The LH process occurs only at surface temperatures 



higher than 440 K for bas al plane sites, and a t even higher temperatures arounc 



step edge and defect sites (jZecho et al 



2mm 



Guttler et al. 



mi 



Zecho et al 



850 K for 



mi). 



The relevance of these results to interstellar dust materials is thus not yet clear. If in- 
terstellar dust materials exhibit these features, it may indicate that in PDRs (where typical 
grain temperatures are between 30-50 K), cher nisorption processes do not contribute signifi - 



cantly to H2 formation on dust grain surfaces (jCazaux fc Tielent 



2004; 



Habart et al. 



20041 ). 



Molecular Hydrogen Formation on Porous Dust Grains 7 
4 EVIDENCE FOR POROSITY IN INTERSTELLAR DUST GRAINS 



The formation and evolution of interstellar dust grains are complex processes which involve 
accretion and chemical reactions of impinging atoms and molecules, grain-grain collision 
and coagulation, photolysis and alteration by U V starlight. X-rays and cosmic rays, erosion 



by sputtering and vaporization (jDraine 



20031). Due to the complexity of these processes. 



there is no complete model that accounts for all the relevant properties of interstellar dust 
grains. The known properties are based on a combination of astronomical observations, 
theoretical studies and analysis of meteorites and solar-system dust. Useful insight about 
the porosity of ii iterstellar dust grains can be obtained from its effect on the observed 



extin c tion curves ( Wrie;htlll987 



I 



.Toneslll988: 



Mathi; 



1996 



Snow fc Wit 



1996 



Foeel fc Leung 



1998; 



Wolff et al 



Wurm fc Blum 



1998; 



20001 : lati et al I 



Schnaiter et al 



lati et al.ll200iri. Interp 



in the upper atmosphere of Earth (jBrownlee 



1999; 



Vadia fc Gupta 



199S; 



Krivova fc IFin 1200 



3 



anetary dust particles and meteorites, collected 



19851 ). were found to be fluffy and porous (Fig. 



P). Although these particles are larger than the typical interstellar dust grains, it is likely 
that similar processes are responsible for their formation. Theoretical and experimental 
simulations of cosmic dust formation, through aggregation of molecular clusters, produce 
fluffy and porous structures with amorphous surface morphology. 

In most calculations of grain temperatures in PDRs, it was assumed that grains are com- 
pact and can be considered as spherical. How ever, several studies considered the tempera- 



tures of non-spherical as well as porous grains (IGreenberg 



1980; 



Voshchinnikov et al 



1999 : iVoshchinnikov et al 



mi 



Blanco Russoletti 



20051 ) . It was found that porous grains 



are significantly colder, in most cases, than predicted by models which use compact spherical 
grains. While models of spherical grains predict grain temperatures in the range of 30 — 50 
K, the actual temperatures of porous interstellar grains may thus be significantly lower. In 
most favorable cases, the grain temperatures may go down to the range of 15 — 30 K for 
large grains and to 10 — 20 K for small grains ([Blanco fc Bussolettilll980| ). 



It was suggested long ago that micro-pores in interstel 



particles and affect the formation of molecules (|Abadi et al 



ar du st grains may trap gas 



19761 ). Recent experiment on 



ice surfaces at low temperatures (10— 25 K) have shown that the porous morp hologv increases 



the desorption energy barriers of adsorbed hydrogen atoms and rnolecules (IPirronello 



1999; 



Manicoe ; 



Viti et al. 



12004 : 



al 



2001 



Roser 



Collings et al. 



et al 



2002 



Hornekaer et al 



2003 . 



200 



3 



Perets et al 



2005 



20041 1. Analysis of some of these experiments shows that 
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on such surfaces would form efficiently only at low surface temperatures of up to 20 



K (IPerets et al 



199fl : 



2005ri. resernb ling the results obtained for olivine and carbon (jKatz et al 



Cazaux fc Tielens 



20041 ). In PDRs, the relevant grain surfaces are expected to consist 



of bare silicate and carbon, not covered by ice mantles. Below we examine the role of grain 
porosity in extending the efficiency range of H2 formation toward higher temperatures. 



5 H2 FORMATION ON POROUS DUST GRAINS 

To evaluate the formation of H2 on porous grains we introduce a suitable rate equation 
model. In this model we make a distinction between the external surface of the grain and 
the internal surfaces within the pores. Atoms adsorbed on the external surface behave as 
in the previous models. They may hop between adjacent sites or desorb. In addition, they 
may enter the pores and slits in the surface. In the internal surfaces, similar diffusion and 
desorption processes apply. However, an atom which desorb from an internal surface, may 
re-adsorb inside the pore, rather than leave the grain. Such an atom may leave the grain 
only after it reaches the external surface, through a series of diffusion and desorption moves 
(Fig. EI). Detailed modeling of these processes requires knowledge on the exact geometry 
of the grains and their surface morphology. However, the most important effects of grain 
porosity can be analyzed by a simple rate equation model. The model, described below, 
assumes high conn ectivitv of the pores , namely, that most of the internal surface of the 



grain is connected (jKimmel et al 



200 ih and is treated as a single surface. The model thus 



consists of two coupled rate equations 



A^ext =F- (^) aN,^, + (^) aN,,, - 2 (^) N^J - WN,., (9a) 



-"ext / ^ '-'in / ^ '-'ext 



iVin = ^ aiVe.t-Pf^ aiVi.-2^ iV,,^ (9&) 

where A'ext and A^in are the populations of H atoms on the external and internal surfaces, 
respectively. The first term in Eq. ()9a|) describes the incoming flux of H atoms. The second 
and third terms in Eq. ()9a|) [and the corresponding first and second terms in Eq. (j9Sj) ] 
describe the diffusion of atoms between the external surface and the pores, where a is the 
hopping rate given by Eq. The numbers of adsorption sites on the external and the 
internal surfaces of the grain are ^cxt and ^in, respectively. In addition, there is a small 
number, S^dge-, of sites, located at the edges of the pores, thus connecting the internal and 
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the external surfaces of the grain. An adsorbed atom moving between the internal and the 
external surfaces must pass through one of these edge sites. The fourth term in Eq. ()9a|l 
corresponds to recombination of atoms on the external surface. Similarly, the third term in 
Eq. (|9&|) corresponds to the recombination of atoms on the internal surface, inside the pores. 
The last term in Eq. ()9a|) describes the desorption of atoms into the gas phase, where W is 
the desorption rate given by Eq. (jlj). 

To examine the effect of porosity on molecular hydrogen formation, we performed nu- 
merical integration of Eqs. under steady state conditions for grains of diameter d = 10~^ 
(cm). The porosity was characterized by the ratio between the internal and external surfaces. 
Sin/ Sext- The formation rate of molecular hydrogen per grain is given by 

i?g.ain = (£^) Ni, + (£) Nl 

In Fig. 121 we show the recombination efficiency of molecular hydrogen on porous grains, vs. 
grain temperature, in the range between 10 — 35 K. Different levels of porosity are used, 
namely 5'in/5'ext = 10, 100 and 1000. The results show that porosity in interstellar grains 
extends the temperature range of high efficiency toward higher temperatures. The extension 
of the temperature range is logarithmic with the level of porosity (Fig. Ej). 

The number of edge sites in the grain is given by S'odge/5'ext = 5 • 10~^. We found that 
although the edge sites control the flow between the external and internal surfaces, the 
recombination efficiency is only weakly dependent on S'edge- In particular, we performed 
simulation with S'edge values 10 times larger and smaller than the value reported above. In 
both cases, the edges of the high-efficiency window were only slightly shifted (not shown). 

In Fig. El we show the number of H atoms on the internal and external surfaces (solid 
and dashed lines, respectively) vs. temperature, for a grain of diameter d = 10^^ (cm) where 
S'in/5'cxt = 100. It is found that for the relevant temperatures there is a large population 
of H atoms in the pores, thus the rate equations are suitable for the calculations presented 
here. The inset in Fig. 0] shows the coverages of H atoms on the internal and on the external 
surfaces. Both coverages are very small. Since the desorption barrier of hydrogen molecules 
on the amorphous carbon surface is lower than the desorption barrier of atoms, the coverage 
of molecules is also expected to be small. This justifies the approximation made in this 
paper, ignoring the population of molecules on the surface and neglecting the LH rejection 
term. 



(10) 
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6 DISCUSSION 

The results presented above show that porosity extends the range of efficient hydrogen 
recombination to higher temperatures. This extension does not reach the grain temperatures 
of 30 - 50 K expected in PDRs and thus cannot, by itself, explain the high abundance of H2 
in these regions. However, porous grains, particularly small ones, are expected to be colder 
and thus may be within the extended temperature range of high efficiency. 

The roughness of interstellar grain surfaces may broaden the distribution of binding en- 
ergies compared to the experimental samples. Thus, some of the adsorbed H atoms may 
be more st rongly bound, further ex tending the range of high efficiency towards higher tem- 
peratures (jCuppen fc Herbstll2005r i. Composite grains, which consi st of both silicates an d 
carbon also exhibit a broader temperature window of high efficiency ( Chang fc Herbsi]l2005 ). 
Grains of diameters smaller than 10~^ (cm) account for most of the surfa ce area of grains in 



interstellar clouds and may thus contribute significantly to H2 formation (jLipshtat fc Biham 



20051). Temperature fluctuations may also reduce the average temperature of small grains 
comapred to large ones. 

In addition to molecular hydrogen, dozens of other molecular species have been ob- 
served in interstellar clouds. Many of these molecules are formed through gas phase pro- 
cesses not involving dust grains. Others are formed on dust-grain surfaces. For example, 
certain organic molecules are believed to be formed on dust grains covered by ice mantles 



(jTielens fc Charnlev 



19971 ). The effect of dust morphology on H2 formation may indicate 



that it plays a role in other gas-grain processes. Moreover, our results may have implications 
on the formation of molecules which have been previously thought to be formed mainly 
by gas phase reactions. For example, recent ob servations of i nterstellar N2 have indicated 



much higher abundances than predicted before (jKnauth 



2004 ). It was suggested that these 



abundances may be explained by gas-grain reactions that had not been considered earlier. 
An analysis of these processes on porous grains may help to understand this new puzzle. 



7 SUMMARY 

We have considered the effect of porosity on the formation of molecular hydrogen on in- 
terstellar dust grains. The analysis was done using a rate equation model, which takes into 
account both the internal and external surfaces of the grain. It was found that porosity 
gives rise to enhanced H2 formation rates. It extends the range of high recombination effi- 
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ciency toward higher temperatures. This may contribute to the understanding of molecular 
hydrogen formation in PDRs, which cannot be explained by models of compact dust grains. 
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Figure 1. Interplanetary dust grain found in the stratosphere, which exhibits the typical porous structure, but at a much 
larger scale. The sizes of interplanetary dust grains range from several nm to 100 fim. Acknowledgement is made to NASA for 
allowing reproduction of this picture from web site |http:/ /stardust.jpl.nasa.gov| 
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Figure 2. Schematic description of H diffusion and desorption in the internal and external surfaces of a porous dust grain. 
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Figure 3. H2 recombination efficiancy for non-porous grains (dashed line) and for porous grains (solid lines) with Sin/Sext = 
10, 100 and 1000. The more porous grains exhibit higher bounds for efficient recombination. 
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Figure 4. The number of H atoms on the internal (soHd Hne) and the external (dashed line) surfaces of a porous grain with 

Sin/Sext = 100 as a function of the grain temperature. Most of the H2 molecules form on the internal surface. The inset shows 
the coverages of H atoms on the internal (solid line) and the external (dashed line) surfaces. These coverages turn out to be 
nearly identical. 



